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Abstract 

Previous work has shown that the macroscopic structure of the theory of quantum 
gravity defined by causal dynamical triangulations (CDT) is compatible with that 
of a de Sitter universe. After emphasizing the strictly nonperturbative nature of 
this semiclassical limit we present a detailed study of the three- volume data, which 
allows us to re-confirm the de Sitter structure, exhibit short- distance discretiza- 
tion effects, and make a first detailed investigation of the presence of higher-order 
curvature terms in the effective action for the scale factor. Technically, we make 
use of a novel way of fixing the total four-volume in the simulations. 
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1 Introduction 



Mundane field-theoretical descriptions of quantum gravity have been undergoing 
a renaissance, thanks to the use of improved renormalization group techniques 
[l] 1 , new ideas about the relation between space and time at short scales [3], and 
a novel way of implementing a lattice regularization of four-dimensional quantum 
gravity, which takes into account both the dynamical and the causal, Lorentzian 
nature of spacetime [4]. 2 New (and partially overlapping [7, 8, 9, 10, 11]) results 
obtained in these approaches give rise to the hope that standard tools from quan- 
tum field theory - adapted to the situation where spacetime itself is dynamical - 
are indeed sufficient to construct a theory of quantum gravity nonperturbatively. 
This is an attractive prospect, since it may imply a large degree of uniqueness, 
with only a small number of parameters needing to be fine-tuned to get to the 
correct theory. 

In this article we will discuss several aspects of the semiclassical limit of the 
research program on "Quantum Gravity from Causal Dynamical Triangulations 
(CDT)", in which a dynamical lattice provides a geometric UV cut-off. In earlier 
papers we have reported that the infrared limit of CDT allows an interpretation as 
the classical solution to Euclidean Einstein gravity with a positive cosmological 
constant [12, 13, 14, 15, 16] (see also [17] for a pedagogical review). Here we 
investigate this limit in more detail, including discussions of discretization effects 
and of the asymmetry between space and time, which appears at the lattice- 
regularized stage of the CDT set-up. 

It is important to bear in mind that what we call the semiclassical limit is a 
truly nonperturbative limit. This means that the tentative continuum limit of the 
path integral is to be found in a region of the bare coupling constant space where 
the entropy of various geometric configurations makes a contribution which is 
as important as the contribution coming from the exponential of the action. In 
lower dimensions, this situation is illustrated by the famous Kosterlitz-Thouless 
transition in the XY model. The XY model is a lattice spin model, whose "spins" 
are two-dimensional vectors of unit length. In two spatial dimensions, this model 
has vortex configurations, with an energy per vortex of approximately 

E = Khi(R/a), (1) 

where k is a coupling constant, R a measure of the linear size of the system and 
a the lattice spacing. Ignoring boundary effects, the centre of the vortex can be 
placed at any one of the (R/a) 2 lattice points. Saturating the path integral (the 

1 building on the original idea of "asymptotic safety" [2] 

2 building on the earlier idea of "dynamical triangulation" as a regularization of quantum 
gravity [5] , but incorporating a concept of micro-causality along the lines of [6] 
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partition function) Z by single-vortex configurations, we obtain 3 



Z = e 



F/k B T 



spin configurations 



e -E[spin]/(fe s T) 



( 




e -[K\n(R/a)]/k B T 



(2) 



We note that the factor (R/a) 2 is entirely entropic, simply arising from counting 
the possible single-vortex configurations, and is independent of any "bare" cou- 
pling constants (the spin coupling k and temperature T). Since the corresponding 
entropy S = ks ln(number of configurations) has the same functional form as the 
vortex energy, we can express the free energy as 



The Kosterlitz-Thouless transition between a low-temperature phase (where vor- 
tices play no role) and a high-temperature phase (where vortices are important) 
occurs when F — 0, i.e. when the entropy factor is comparable to the Boltz- 
mann weight of the classical energy. At this point we are far from the naive 
weak coupling limit of the lattice spin theory, which is just a Gaussian free 
field. Instead, the continuum field theory associated with the transition is the 
sine-Gordon field theory at the coupling constant value where it changes from a 
super-renormalizable theory to a renormalizable one. 

The situation in four-dimensional CDT is analogous, in the sense that the 
nontrivial structure of the phase diagram reported in [19] results from a compe- 
tition between the entropy of configurations and the action, precisely as in (3). 
The analogy goes even further: when written as an effective action for the global 
scale factor (which plays the role of an order parameter for the gravity case), 
and in the region of phase space we have identified as possessing a meaningful 
classical limit, the free energy has the same functional form as the classical ac- 
tion, but with the opposite sign (corresponding to F < in (3)). This makes it 
important to understand the semiclassical nature of this phase, and its relation 
to continuum physics. The investigation of this issue is the main purpose of the 
present article. 

The remainder of this article is organized as follows. In Sec. 2 we recall the set- 
up of causal dynamical triangulations and introduce a new way of (approximately) 
fixing the total four-volume of spacetime in the simulations. Using this new 
prescription, we re-analyze three-volume distributions in Sec. 3 and examine the 
probability distributions for different values of discrete three- volumes. The latter 
allows us to exhibit and quantify lattice artifacts for small three-volumes, which 

3 Our present discussion is merely qualitative and meant to highlight the competition between 
entropy and Boltzmann weights; exact treatments of the Kosterlitz-Thouless transition are given 
in many textbooks, see, e.g. [18]. 



F = E-ST = (k- 2k B T) \n(R/a). 



(3) 
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occur close to the beginning and end of our quantum universe. In Sec. 4 we 
introduce a refinement of the spatial slicing, associated with connected layers of 
the different types of simplicial building blocks. This allows us to analyze the 
corresponding volume distributions separately and compare their relative scaling 
behaviour as a function of the bare coupling constants. In Sec. 5 we return to 
the task of reconstructing the effective action for the three-volume fluctuations 
from the computer measurements. In a nontrivial extension of earlier work we 
look for evidence for the presence of corrections to the effective action associated 
with curvature-squared terms in the continuum. We summarize and discuss our 
findings in Sec. 6. The appendix contains a extension of the analysis of Sec. 5 to 
the situation with a refined spatial slicing. 

2 Causal Dynamical Triangulations 

In this section, we will review briefly some key ingredients of our approach and 
give an outline of the methods used. The basic motivation behind quantizing 
gravity via CDT was explained in [20], and the current set-up described earlier 
in [4, 13, 16], to which we also refer for further technical details. 4 

Assuming that a path integral representation of quantum gravity exists, the 
basic idea is to provide it with an ultraviolet cut-off by using piecewise linear 
("triangulated") geometries in the quantum superposition. Assigning lengths to 
the one-dimensional edges of such a triangulation fixes its geometry completely, 
without the need to introduce coordinates [22] , thus also avoiding the redundan- 
cies of the usual continuum description of curved geometries. For our purposes, 
we further restrict the piecewise linear manifolds to those which can be obtained 
by gluing together two specific types of building blocks ( "four-simplices" ) with 
prescribed edge lengths. The typical edge length serves as a geometric UV cut-off 
and is a measure of the fine-grainedness of the geometry. 

Unlike its Euclidean counterpart, CDT also employs a proper-time foliation, 
with respect to which the topology of space is not allowed to change. Admissi- 
ble geometries are those which can be constructed by first triangulating spatial 
slices of constant proper time, which for simplicity we assume to have the topol- 
ogy of S 3 . Each three-slice is assembled from identical building blocks, namely, 
equilateral tetrahedra whose spacelike edges all have length a s (the spatial lattice 
spacing or UV cut-off). The gluing and the number of tetrahedra in each slice 
are arbitrary except for the overall topology and the imposition of local manifold 
constraints. 5 The next step consists in connecting neighbouring spatial slices by 

4 An up-to-date reference describing CDT and key results obtained in it is [21]. 
5 Note that neither of these constraints on individual path integral histories will necessarily 
survive in the continuum limit. Previous results from CDT demonstrate that it can and does 
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timelike edges of length a t (the lattice spacing in time direction), in a way con- 
sistent with the "filling" of the spacetime sandwiches by (up to rotations) two 
types of flat, Lorentzian four-simplices, to wit, 

• simplices of type (4,1), which have four vertices in common with the spa- 
tial slice at (integer proper) time n (thus spanning one of the equilateral 
tetrahedra making up the slice) and one vertex with the slice at time n+1. 
Time-reversed building blocks which share one vertex with slice n and four 
vertices with slice n+1 are "of type (1,4)". 

• simplices of type (3,2), which have three vertices in common with the spatial 
slice at (integer proper) time n (thus spanning one of the equilateral trian- 
gles contained in the slice) and two neighbouring vertices with the slice at 
time n + 1, spanning a spacelike edge there. Time- reversed building blocks 
which share two vertices with slice n and three vertices with slice n + 1 are 
"of type (2,3)". 

These four-simplices are glued pairwise along their three-dimensional "faces", 
forming a layered, four-dimensional simplicial manifold of topology S 3 x [0, 1] in 
the manner just described. In such a piecewise linear geometry the curvature 
is concentrated at the two-dimensional subsimplices (the triangles) and induces 
nontrivial rotation angles on vectors which are parallel-transported around them. 

Two neighbouring spatial slices labeled by integers n and n + 1 are separated 
by a proper-time distance a t in the sense that each timelike edge connecting 
the two hypersurfaces has this length. Instead of using two lattice spacings a s 
and a t , we usually work with a single a := a s > and the dimensionless ratio 
a := —a 2 /a 2 . Expressing a in terms of squared lengths allows us (i) to start out in 
Lorentzian signature, where af < 0, and thus a > 0, and (ii) to perform for each 
triangulation an analytic continuation of a in the lower-half complex a-plane to 
real, negative a, resulting in a piecewise linear geometry of Euclidean signature 
and length assignments af = \a\a 2 to the edges that in Lorentzian signature used 
to be timelike. This "Wick rotation" relies on the foliation in proper time. It 
has no obvious correspondence in the metric continuum formulation (see [23] for 
a discussion). Under this map, the Einstein-Hilbert action of a given Lorentzian 
piecewise linear geometry, the so-called Regge action Sl [22], changes according 
to Sl(cx) — > iSE^—a), where Se(— «) is the Regge action for the corresponding 
piecewise linear Euclidean geometry with the length assignments a 2 = \a\a 2 . A 
necessary condition to ensure that the Euclidean four-simplices are nondegenerate 

happen that the system is driven dynamically to a quantum configuration which is no longer 
of topology S 3 x [0, 1] and no longer resembles a four-dimensional manifold. This reflects the 
nontrivial interplay of "energy" and "entropy" already alluded to in the previous section, which 
the nonperturbative CDT formulation allows us to capture. 
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is \a\ > 7/12, which will be assumed in what follows. Once the rotation to 
Euclidean signature has been performed, we redefine a — > —a, Se(-oc) — > Se{&) 
for simplicity of notation. 

In our construction of the path integral, we start out with the set of causal 
piecewise linear Lorentzian geometries described above. In order to perform the 
sum over these histories, we then rotate each of them to Euclidean signature, so 
that it will contribute with a real weight exp(— Se)- More precisely, we have 

^ = -(/ t0 + 6A)iV + / t4 (iVi 4 ' 1) + ^ (4) 

where N , ivf ' 1} and 7V"j 3 ' 2) are the total numbers of vertices (zero-simplices) , of 
four-simplices of both type (4,1) and type (1,4), and of four-simplices of types 
(3,2) and (2,3) respectively. The Regge action [22] takes the particularly simple 
form exhibited in the bottom line of eq. (4) because each CDT contains only 
two geometrically distinct types of building blocks. The parameter k,q in (4) is 
proportional to the inverse bare gravitational coupling constant, is a linear 
function of the bare cosmological and the inverse bare gravitational constant, 
while A is an asymmetry parameter, which for given re 4 and k can be related 
to a. It is normalized such that a = 1 (the case of equilateral Euclidean four- 
simplices) corresponds to A = 0. 

For technical reasons we perform the numerical simulations at fixed four- 
volume, which in practice is usually realized by fixing N^' 1 ^ to some target volume 
V4. The geometry is updated by using local Monte Carlo moves (c.f. [4, 13, 
16]), which in general will change the four- volume. In past papers this was 
taken care of by modifying the action by a linear term according to Se Se + 
e|(A^ 4 '^ — V4) j , with a parameter e controlling the range of fluctuations, and 
measurements collected only when JV^ 4 ' 1 ^ = V4 exactly. Fixing JV^ 4 ' 1 ^ in this 
delta function-like manner has the side effect of generating a zero mode in the 
correlation matrix {SV3(n)SV 3 (n')) of fluctuations in the spatial t/iree-volume V 3 
in time, an important quantity in investigating the dynamical behaviour of the 
quantum universe [16]. Projecting out the zero mode, which is necessary for 
inverting this matrix, leads to an inconvenient mixing of the remaining modes. 
To get around this in the present paper, we no longer require that 
exactly in measurements, but instead let it fluctuate according to the modified 
action Se + ^iV"] 4 ' 1 "* — V4) 2 , which eliminates the zero mode. The advantage of 
using a quadratic term is that it can be treated easily alongside other terms in 
the action. 

In a given simulation, the parameter k 4 of the model is fixed by requiring that 
(N^' 1 ^) = V4. The data presented below were taken at e = 1CT 5 , but we checked 
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that the results are essentially unaltered for e = 2 x 10 5 and e = 5 x 10 6 . 
Furthermore, we also checked (at e = 2 x 10~ 5 ) that using the total number N± of 
four-simplices instead of N^ 4 ' 1 ^ as the target volume does not change the situation 
either. Ideally one would like to make e as small as possible, but this must be 
balanced against the fact that smaller e- values increase the auto-correlation time. 

Depending on the values of the bare couplings in (4), the CDT theory will 
appear in one out of three different phases A, B or C (see Fig. 6 for a depiction 
of the phase diagram in the /t -A-plane, and [19] for a detailed description). In 
the present paper we will concentrate on phase C where an extended de Sitter 
universe has been observed [15, 16, 19]. In this phase we can construct an effec- 
tive semiclassical action for the scale factor of the universe. Changing the bare 
coupling constants n and A will affect the effective coupling constants appearing 
in this action. 

3 Volume distributions reloaded 

An interesting observable investigated previously is the typical "shape" or "vol- 
ume profile" of the quantum universe, more precisely, the average distribution 

(N(n)) := (N^\n)) (5) 

of spatial volume as function of the discrete proper time n. 6 What we found 
is that individual configurations in the well-behaved "de Sitter phase" have a 
characteristic shape, consisting of a (spatially extended) bulk contribution or 
"blob" and a separate, thin "stalk", whose spatial volume for all times stays close 
to the prescribed minimal cut-off size. Fig. 1 shows an example of the volume 
distribution N(n) of a typical, individual path integral configuration, compared 
to the average distribution (N(n)}. 

If we compactify the time direction, as is usually done in simulations, the 
action is symmetric under time translations, permitting us to move the position 
of the "centre of volume" along the time direction. To determine the average 
distribution we shift this position for every individual configuration such that its 
centre coincides with a common reference point n on the time axis. For each 
shifted configuration we measure the three- volume distribution N(n) obtained in 
numerical simulations at some fixed V4 as a function of the shifted time variable 
n, which can be averaged meaningfully over many independent configurations. In 
[14] it was shown that inside the blob region the average distribution obeys the 
characteristic scaling relation 

(N(n)}=A s Vt 1/dH cos dH - 1 (r), (6) 

6 Equivalently, one can work with the number N 3 (n) of three-simplices, where N^ ,1 \n) = 
2N 3 (n) for all interior constant-time slices of the triangulation. 
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Figure 1: Example of a volume distribution N(n) for a specific path-integral 
configuration (red), compared to the average distribution (N(n)) (blue). On the 
scale of the plot, the latter cannot be distinguished from a fit to the theoretical 
curve of relation (6). 
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where r = A t n/V 4 H represents proper time and dn = 4 (within numerical ac- 
curacy), and A t and A s are numerical constants. In the stalk part the average 
volume is practically constant and independent of V4. The fact that the numer- 
ical data reproduce the distribution (6) with great accuracy remains practically 
unchanged when the measurements are performed with the new method. 

The interpretation of the above results is that the quantum system generates 
the three-volume distribution of a round four-sphere (i.e. Euclidean de Sitter 
space) with four- volume V4, with time steps proportional to the (cosmo logical) 
proper time separating the spatial hypersurfaces. The stalk is present only be- 
cause (i) we have fixed the time period T to be larger than the proper-time 
extension of the universes created, and (ii) we enforce a minimal three-volume 
of 5 tetrahedra per spatial slice, corresponding to iV min (n) = 10, in keeping with 
the simplicial manifold character of the spacetime. Although it appears from the 
simulations that the dynamics of the system wants to drive the three-volume to 
zero at the two ends of the blob, we prefer to maintain these kinematical restric- 
tions (and thus the stalk region), because they allow us to monitor fluctuations 
in the time extent of the blob, which would be suppressed if we tried to adjust 
the time interval T to match the blob exactly. 

The measurements presented below were taken for the coupling constant val- 
ues (k , A) = (2.2,0.6), four-volume V 4 = 160 000 and time period T = 80, and 
exhibit a behaviour typical for systems inside the de Sitter phase. The reference 
time has been fixed to no = 40.5 and the distribution symmetrized with respect 
to n — > 81 — n. Comparing the volume distributions with relation (6), the range 
of the blob is roughly 17 < n < 64 in terms of integer times, whereas the stalk 
region is located in 1 < n < 17 and 64 < n < 80. 

In previous work [14, 15, 16] we have matched the average (N(n)) to the clas- 
sical solution N c i(n) of a (discretized) mini-superspace action. This semiclassical 
interpretation is corroborated further by new measurements made of the prob- 
ability distributions V n (N), which for a given time n describe the probability 
that the three-volume N(n) is equal to N. If the semiclassical picture is indeed 
correct, these distributions should be approximately Gaussian with a mean equal 
to (N(n)) and a dispersion a 2 , = (N(n) 2 ) — (N(n)) 2 . This is exactly what we 
observe in the range 25 < n < 56, well inside the blob, as shown in Fig. 2 (top). 

Inside the stalk, the situation is completely different. Not only the average 
(jV(ra)), but also the three- volume distribution V n (N) is unchanged when we 
move along the stalk. In addition, we observe short- distance lattice artifacts, 
in the form of a split of the distribution of three-volumes into three separate 
families. This is clearly visible in Fig. 2 (bottom), which shows the average 
three- volume distribution in the stalk for 1 < n < 17. Inside each family the 
discrete three-volumes differ by 6, such that we have three sets {10, 16, 22, . . . }, 
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Figure 2: Top: probability distribution V n (N) of three-volumes at fixed time 
n = 29; bottom: probability distribution V n (N) in the stalk (for every n < 17); 
the data fall into three families (colour-coded in the graph), each one with a 
different behaviour. 
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{12, 18, 24, ... } and {14, 20, 26, . . . }. 7 Although the presence of the stalk - as 
we have argued above - is merely a consequence of our kinematical set-up, it 
cannot simply be ignored, but has to be taken into account in the calculations, 
for example, when analyzing the covariance matrix of three-volume fluctuations 
and its inverse. In order to distinguish short-distance artifacts from short-distance 
physics, it is therefore important to understand the behaviour of the transition 
region between stalk and blob as best possible. 

An analogous split into three sets of the three-volume distribution for small 
volumes is also observed in the region 17 < n < 25, which smoothly joins the 
stalk and which we have classified as "inside the blob", but where (N(n)) is still 
relatively small. We find that for 10 < N(n) < 200 the distribution again splits 
into three families. An example for n = 22 is shown in Fig. 3 (top); the situation 
is similar for 18 < n < 24. If in the same transition region one considers higher 
values of N(n) > 200, the split between the three families disappears, but one 
also notes that the distributions are highly asymmetric and non-Gaussian, as 
again illustrated for n = 22, depicted in Fig. 3 (bottom). 

We conclude that lattice artifacts are visible for spatial slices of sizes up to 
A" 3 = 100 tetrahedra (corresponding to N = 200). This may sound like much, but 
in reality corresponds to rather small linear distances. This can be understood 
by comparing the total volume of 100 spatial tetrahedra, 100 x ^| a 3 ., with that 
of a regular three-sphere, 2n 2 R 3 . If we arranged the tetrahedra to approximate 
such a sphere, we would have R ~ a s , i.e. only a single lattice spacing, with the 
antipodal distance on the three-sphere being roughly equal to 3a s . Viewed like 
this, it appears rather encouraging that the observed lattice artifacts vanish so 
quickly as function of the three-volume. 

Quantifying the short- distance artifacts for one quantity gives us a good idea 
of the scale at which they occur for a given four-volume, but still leaves us with 
the difficulty of disentangling them from genuine quantum-gravity signatures in 
this and other observables. For example, the non-Gaussian character of the dis- 
tribution of fluctuations around the sphere (6) at small A 3 (n) observed above 
could indicate new short-distance physics, say, along the lines suggested in the 
asymptotic safety scenario [24]. We presently do not know how to relate the 
deviation from the four-sphere at small scale factor a(r) described there to our 
explicit small-A(n) "observations", but it would clearly be interesting to do so. 

4 Refining the spatial slicing 

An alternative description of the four-dimensional simplicial manifolds contribut- 
ing to the sum over histories is given in terms of their dual structures, where we 

7 Note that by definition N(n) is always an even number. 
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Probability distribution for the transition region 
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Probability distribution for the transition region 
0.009 , , , , , 




Figure 3: The probability distribution P22{N) from the transition region near the 
end of the blob. For small N the distribution splits into 3 families (top). For 
large N the split disappears, but the distribution is highly asymmetric (bottom, 
no colour-coding for the three families). 
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n 



Figure 4: The averaged distributions (N(n)) and (^( n + |)) combine into a 
single curve after performing a suitable relative rescaling. 

assign vertices to the centres of four-simplices and links to the three-dimensional 
faces between adjacent four-simplices. In the absence of boundaries of the mani- 
fold, each dual vertex will have five neighbouring ones, connected by dual links. 
This dual picture suggests a further subdivision of the original (integer) time 
steps of the discrete triangulation, as we will now explain. 

Integer steps n provide natural units for timelike paths in the original lattice 
if the paths are taken to run only along the timelike 8 edges of the four-simplices, 
which form one-dimensional 'connectors' between pairs of adjacent constant-time 
layers made entirely from spatial tetrahedra. Constructing analogous paths on 
the dual lattice, at least four steps are required to connect a vertex dual to a 
(4,l)-simplex at time n with a vertex dual to a (4,l)-simplex at the next time 
step n+l. Labelling the dual vertices by their associated simplex types, such a 
path is given by a sequence (4, 1) — > (3, 2) — > (2, 3) —> (1, 4) — > (4, 1). In addition, 
we have observed that all dual vertices of types (3,2) and (2,3) between times n 
and Ti + l and the dual links connecting them form a single closed, connected 
graph. This makes it natural to assign a time n + 1/2 to this layer and study the 
properties of the distribution of 

N(n + 1/2) :=Af' 2 V + l/2), (7) 
8 "timelike" refers to the character of the geometric elements before the Wick rotation 
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n 



Figure 5: Combining the averaged volume distributions (N (n)}, (iV (n + |)) and 
(N (n+ I)) into a single one. 

where N± ' (n + 1/2) is the total number of (3,2)- and (2,3)-simplices located 
between discrete times n and n + 1. It is not surprising that this distribution is 
similar to (N(n)), in fact, by a simple rescaling 

N(n + 1/2) = P N(n + 1/2) (8) 

(where the constant p will depend on the values of the couplings), one can achieve 
that in the blob range the combined distribution (N(n)), with n now running over 
both integer and half-integer times, is well approximated by a single, smooth curve 
according to formula (6). Note that this unification is not global, in the sense that 
a different rescaling is needed in the stalk part, which is not surprising in view of 
the different dynamics in this region. Fig. 4 shows the resulting distribution for 
the standard choice of bare couplings, Kq = 2.2, A = 0.6. 

One can push this line of argument further by performing an even finer dis- 
crete subdivision of time. As argued above, going forward in time along dual 
links necessarily takes one through a sequence of different types of four-simplices. 
Nothing new can be learned by distinguishing between a layer of (4,l)-simplices 
and that of the adjacent (l,4)-simplices, because they are in one-to-one corre- 
spondence (each interior spatial tetrahedron at integer n is shared by exactly one 

(4.1) - and one (l,4)-simplex). By contrast, there is no such relation between the 

(3.2) - and the (2,3)-simplices. This has motivated us to divide each integer time 
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step into three and study volume distributions separately for the sets of (4,1)-, 
(3,2)- and (2,3)-simplices contained in each 'sandwich' [n,n+ 1]. One finds again 
that they can be combined into a single, universal distribution by defining 

N(n+ 1/3) := 2pN {3 ' 2 \n + 1/3), N(n + 2/3) := 2piV (2 ' 3) (n + 2/3), (9) 

with the same p as used before. In (9), N ( - 3 ' 2 \n + 1/3) and N^^n + 2/3) count 
the number of (3,2)- and (2,3)-simplices, which collectively have been assigned 
the time labels n + 1/3 and n + 2/3, respectively Fig. 5 illustrates that the 
sphere fit continues to work beautifully also with respect to this time subdivi- 
sion. These results underscore that the detailed choices we make at the level of 
the individual building blocks, in the present example the identification of mi- 
croscopic time steps, bear little direct relation to the macroscopic aspects of the 
semiclassical emergent geometry. This may have been anticipated because of the 
nonperturbative nature of this limit in CDT quantum gravity. 

Our next step will be to repeat the above measurement for different values 
(/to, A) of the bare coupling constants, and to study systematically how the scaling 
parameter p behaves as a function p = p(no, A). We are particularly interested in 
the behaviour inside the de Sitter phase (labelled C in Fig. 6, top) as we approach 
one of the critical lines, which separate phase C from phases A and B. The set of 
points in the coupling constant-plane where measurements were taken lies along 
the T-shape inside phase C. 

As is apparent from Fig. 6, bottom right, there is practically no dependence 
of p on A for fixed Kq, despite the fact that the quality of the fits to the functional 
form (6) deteriorates for A < 0.3. On the other hand, leaving A fixed and varying 
kq, the dependence becomes quite pronounced, with p increasing monotonically 
as the A-C phase transition is approached. However, with the current data it is 
not really possible to nail down uniquely the functional form of the dependence, 
and thus the nature of the phase transition. Let us illustrate this point by drawing 
two different three-parameter fits through the data points of Fig. 6, bottom left. 
(Note that the error bars of the data points themselves are too small to be visible 
in the figure.) For fit 1 (red curve) we choose the form 

p = A(k c - k )- c , (10) 

for which we have determined the free parameters as 

A = 1.67 ± 0.03, /eg = 4.76 ±0.02, C = 0.82 ± 0.03 (x 2 « 1). (11) 

Given the excellent agreement with the data and the small error bars, one is 
tempted to conclude that the transition is characterized by a fractional critical 
exponent, a result which usually indicates a higher-order phase transition. How- 
ever, such a conclusion would be premature. The error bars of such fits can be 
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0.5 1 1.5 2 2.5 3 3.5 4 



Figure 6: Systematic study of the dependence of the parameter p on the couplings 
Kq (bottom left) and A (bottom right); the points in the phase diagram at which 
the measurements were made lie along the T-shaped line in phase C (top figure). 



misleading, since they refer to a specific functional form, which is unlikely to 
be correct when we move away from the critical line. As an alternative, let us 
assume that 1/p is analytic in and goes to zero linearly as we approach the 
critical value Kq, instead of with a fractional power. Matching the number of free 
parameters, the corresponding fit 2 (green curve) has the form 

1/p(ko)=A + B*Kq + C*i4, (12) 

and its parameters have been determined as 

A = 2.10 ± 0.02, B = -0.30 ±0.02, C = -0.026 ± 0.003 (x 2 « 2), (13) 

corresponding to a critical value Kq = 4.86 ± 0.22. We note that this fit is almost 
as good as the previous one (with \ 2 ~ 2 rather than x 2 ~ 1), and the two curves 
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are barely distinguishable on the plot. However, it leads to the rather different 
conclusion that the critical exponent is —1 rather than —0.82, indicative of a 
first-order transition. The fact that the exponent —1 lies more than 5 standard 
deviations away from —0.82 reiterates our earlier assertion that the error bars 
given with the individual fits cannot be taken seriously as the only source of 
error. Deciding between these two possibilities will require taking data closer to 
the critical line, which is quite difficult from a simulation-technical point of view. 

To summarize, the data presented above do not allow us to distinguish be- 
tween a first- and a higher-order transition. The possibility of a higher-order 
A-C transition is intriguing, since it is exactly by changing k that we encounter 
a transition scenario a la Kosterlitz-Thouless as described in the introduction, 
however, it is not supported by the results reported in [19]. There we observed 
an abrupt transition and clear evidence of a hysteresis when changing the cou- 
pling constants, which we interpreted as good evidence of a first-order transition. 
The results reported above do not really contradict this hypothesis, although they 
leave the door ajar to the possibility of a higher-order transition. More work is 
needed to settle this question decisively. 

5 Effective action: curvature corrections 

As mentioned above, the measured averaged volume distribution (N(n)) closely 
matches the behaviour of the same quantity derived from a discretized mini- 
superspace action. In [16], we made an effort to push this line of argument further 
by trying to determine the form of the effective action from measurements of 
the covariance matrix of the three- volume fluctuations 5N(n) = N(n) — (N(n)). 
Extracting these data was more subtle, in part due to the presence of a zero-mode, 
whose origin was the constraint of keeping the number of (4,l)-simplices constant. 
This mode had to be projected out before one could invert the fluctuation matrix 
and use it to reconstruct an action. An unpleasant feature of this projection is 
its mixing of effects from the stalk with "bulk physics". One new ingredient in 
the current work is to lift the constraint of a constant number of (4,l)-simplices, 
and thus avoid the zero-mode problem. 

In this new setting, we study the covariance matrix C n>n > = (8N(n)5N(n f )) 
and invert it to a matrix S ntTl r which in principle describes the (nonlocal) effective 
action resulting from fluctuations of the three-volume around the semiclassical 
solution N c i(n) = (N(n)) (see [15] for details). Because we have added a term 
e (J2 n N(n) — V4) 2 to the original Einstein-Hilbert action, each matrix element of 
S n>n > receives a contribution 2e. After subtracting this shift, we find that all matrix 
elements outside the diagonal and the neighbouring sub- and superdiagonals are 
zero up to numerical noise. To quadratic approximation in the fluctuations we 
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have 



d 2 s eff 



Sn ' n ' dN(n)dN(n') 



(14) 

N(n)=N cl (n) 



Our observations then suggest that the effective action is quasi-local in time and 
can be expressed as 

Se ff = ^ (f(N(n),N(n + l)) + V(N(n)j) + e( £ N(n) - V,)\ (15) 

n n 

where the function / and the potential V need to be determined from the data 
and T is an overall constant. The simplest function compatible with the observed 
scaling is 

><*■»> = l^F- < 16 > 

Assume now that the three- volume (equivalently, the third power a 3 (n) of the 
scale factor of the universe) behaves according to 

{ N(n)) = A,V: / 'H<^) (17) 

for some function H, where no denotes the location of an origin chosen along the 
discrete time axis. Comparing with (6), this scaling is consistent with our data 
and we have 

N d (n) = A S V? /4 cos 3 (A(n - n )/V 4 1/4 ), A s = ^A t , (18) 

where the last equality is required by normalization. Converting the function 
f(N(n + 1), N(n)) found in (16) to a continuum expression yields 

/Wn + W) ^(^) 2 , (19) 

where t oc na t /V^ 4 is a dimensionful time variable, and a t the lattice spacing in 
time direction. 

We recognize (19) as the kinetic term of a minisuperspace action for a spatially 
homogeneous and isotropic universe. One could consider a variety of corrections 
to this term. More specifically, we are looking for corrections related to the short- 
distance behaviour of the theory. Taking guidance from the continuum theory of 
a homogeneous, isotropic universe with (Euclideanized) four-metric 

ds 2 = dt 2 + a{t) 2 dtt 2 {3) , dtt 2 {3) = d9 2 + sin 2 Q{d$\ + sin 2 fa d<f) 2 2 ), (20) 
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the most general minisuperspace action containing both the Ricci scalar and 
a curvature-squared term 9 (the latter with coupling constant u), as well as a 
cosmological-constant term is of the form 




( 



a a 




(21) 



We note the presence of fourth-order time derivatives, whereas the higher-order 
spatial derivatives in such a universe are converted to inverse powers of the scale 
factor a(t). The i?-term multiplied by the inverse Newton constant in the inte- 
grand of (21) contains also a potential part a oc N 1 ^ 3 , which in earlier work [16] 
was matched successfully to computer measurements. The fact that correspond- 
ing terms in the effective action (15) appear with the opposite sign compared to 
(20) has to do with nonperturbative contributions arising from integrating over 
the remaining degrees of freedom in the path integral [13, 14, 16, 21], and is 
therefore another feature of the nonperturbative nature of CDT's semiclassical 
limit in phase C, 

For the purposes of our present investigation, do we observe any trace of the 
terms contributing to the i? 2 -term in our effective action? With regard to the 
purely diagonal potential term V(N(n)) in (15), we would expect higher-order 
corrections 10 in powers of iV _2//3 , that is, 



V(N(n)) = -X eff N{n) + fiN{n) 1/3 + ^N{n)- 1/3 + (^(n)- 1 ) (22) 



to the order we are considering, where the first term is a Lagrange multiplier 
term. 11 Next, let us turn to the kinetic term in S e g. Since we have already 
seen that only nearest-neighbour terms contribute appreciably to the measured 
covariance matrix, we can immediately conclude that a term corresponding to 
a(t)a(t) 2 - whose natural discretization would contain a term like (N(n + 1) — 
2N(n) + N(n — 1)) 2 /N(n), and therefore next-to-nearest-neighbour interactions 
- is absent. This is our first indication that in the region of phase space under 
investigation, there is no appreciable contribution from a squared-curvature term 
matching the one in the continuum expression (21), when re-expressed as function 
of the three-volume N. 

Let us look at another kinetic term expected to contribute to the action at 
this order, a 4 /a oc N 4 /N 3 , to illustrate how a more quantitative analysis needs 

9 For a classical universe with metric (20), the a priori distinct terms R 2 , R^R^ and 
R i iv P <yR^ uf " y in the action are all proportional. 

10 Such a power expansion should be valid as long as N(n) does not become small, i.e. suffi- 
ciently far away from the "beginning" and "end" of the universe. 

11 From expression (22) it is clear that if the first two terms should both contribute to leading 
order in the continuum limit, at least one of the multiplying constants needs to scale nontrivially 
as a function of N. This is borne out by relation (25) below. 
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to proceed and what the potential implications are for the continuum limit of the 
theory, if we manage to construct one. Pretending for the time being that this 
term is the only correction to the finite-difference expression (16) at this order 
leads to a modified ansatz 

/(!t ,,)=<£^(l + 6 (£ZVy + ...Y (23) 
x + y y \x + yj J 

As mentioned earlier, the higher-derivative expansion takes the form of a power 
series in TV -2 / 3 or, equivalently, in inverse powers of the square root \fV± of the 
discrete four- volume (the number of four-simplices). It is our task to determine 
from the data whether or not such subleading terms are present in the effective 
action. However, let us emphasize that the presence of higher- derivative correc- 
tions at the regularized level, like those associated with the coefficients £1 and 
£2 in (22) and (23), does not imply that such terms necessarily survive in the 
continuum limit. For example, if in the simulations we observed a coefficient £2 
which was independent of the four- volume V4 of the universe, we could convert the 
dependence on 1/ y/V^ of the corresponding term in the action into a dependence 
on the lattice spacings a t and a s in the time and spatial directions according to 

V 4 cont oc a t a 3 s V 4 , (24) 

or a 4 ^ in shorthand notation 12 . In this case the term (x — y) 2 /(x + y) 2 in (23) 
would be proportional to a 2 and simply drop out in a standard scaling limit where 
we keep Vf ont constant while taking V4 to infinity and the lattice spacings a t and 
a s to zero. Likewise, if we wanted to add such a higher-derivative term by hand 
already in the bare action, such that it survived in the classical limit, we would 
need to assume an appropriate nontrivial lattice dependence of the coefficient £ 2 
in front of the discretized higher-derivative term. However, we will not follow this 
latter route in the present work. This does not imply that no higher-curvature 
terms may appear in the effective action, it only means that no tunable coupling 
constant is associated with such terms, and that the coefficient £ 2 is determined 
purely from the entropy of microstates. 

In the region where we can presently perform reliable measurements, £2 does 
not display any significant scale dependence. However, it could in principle pick 
up such a dependence when we move closer to a potential ultraviolet fixed point 
by changing the bare coupling constants, along the lines described in [16]. In this 
way, in the vicinity of a nontrivial ultraviolet fixed point (in principle infinitely 
many) higher- derivative terms can play a role. In view of this situation our task 

12 The notation "a" for a generic lattice spacing should not be confused with the same notation 
for the scale factor a(t). 



20 



is to first identify potential higher-derivative terms in the effective action, and 
then study the scaling behaviour of the associated coupling constants. 

Let us determine the effect of the presence of a ^-dependent term on the 
fluctuation matrix S njTl i. We can use the fact that our tentative discretized action 
should reproduce the "observed" N d (n) of (18) as an extremum. Extremality is 
satisfied if 

Q ( A 2 fX\ 2 / 3 \ 

V(N(n)) = - l-^=N(n) + ^-J A^N^n) + 0{N-^{n))\ . (25) 

To this order the result does not depend on the parameters £1 and £2- Note that 
A t can be fully determined from the data for N c i(n). Furthermore, the additional 
term proportional to e gives no contribution to the classical solution, since the 
bare cosmo logical constant k 4 was chosen to achieve exactly that, as explained 
at the end of Sec. 2. 

We can now use (15) to derive a prediction for the inverse covariance matrix, 
namely, the matrix S n y of second derivatives at N(n) = N c i(n) (see (14)), pro- 
vided that the fluctuations around N d {n) can be approximated by working up 
to second order in the 5N(n) (see [16] for a detailed discussion). It is a simple 
exercise to show that the matrix elements (after again eliminating a constant 
2e-shift) satisfy 

N d (n + 1) 1 N d (n-1) _ 1 , 

S n , n + Nd{n) S n , n+1 + jVd(n) Sn.n-1 " f V (n) (26) 

where 

/q\ 2/3 

V"(n) = V"(N d (n)) = - ^-J A 8 t /3 N d 5/ \n) + 0(N'^(n)). (27) 

As long as we keep only the leading term in (27), the right-hand side of (26) is 
independent of £ 2 - However, the elements S njTl and ^^±1 can (and do) depend 
on £2 at this order, provided these contributions cancel on the left-hand side of 
(26). This is illustrated by the expression for S n;7l ±i, which after subtraction of a 
2e-contribution is given by 

1 8N d (n)N d (n±l) ( ( N d {n ± 1) - N d {n) \ 2 \ 

n ' n±1 ~ T (N d (n) + N d (n ± !))» + ^ { Ncl (n ± 1) + N d (n) J + ' ' ) ^ 

In Fig. 7 we compare this prediction with the data for the standard choice 
(/■Co, A) = (2.2, 0.6) of bare constants. Notice that the term proportional to £ 2 is 
0{V^ 1 ^ 2 ) with respect to the leading term, exactly as in the "potential" V"(n). 
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Figure 7: Determining the parameter £2 for the standard choice of couplings. 
Top: finding the optimal value of £2 by minimizing the error of 1/r in the fit for 
the subdiagonal S n , n +i (best fit denoted by the green square). Bottom: the best 
fit in a log-log plot of \S n ^ n+ i — 2e\ versus the averaged three- volume {N c i(n) + 
N d (n + l))/2. 
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Fig. 7 shows the results of the analysis including the leading correction to the 
kinetic term. The best fit is obtained for £ 2 — 0.327 ± 0.010. This concludes 
our discussion of how to search for and analyze higher-order correction terms 
to the effective action. For such terms to have a continuum interpretation as 
a curvature-squared contribution like in the continuum action (21), we need to 
show that a whole bunch of coupling constants like £2 exhibit the correct scaling 
behaviour as we follow their flow when we approach the critical phase transition 
lines. Our present Monte Carlo algorithm still needs improvement to deal with 
the critical slowdown appearing close to the B-C transition line, which in this 
context has our primary interest [16, 10]. 

In the appendix we repeat the above analysis for the system with subdivided 
time steps, analogous to what we did in Sec. 4 for the volume distribution. This 
turns out to be self-consistent, in the sense that integrating out the substructure 
at the discretized level produces a result compatible with what was derived for 
integer time steps in the present section. 

6 Summary and discussion 

In the work presented above, we have analyzed several properties of the semiclas- 
sical limit of CDT quantum gravity, a limit which is inherently nonperturbative 
in nature. We used a modified simulation method which replaces the fixing of the 
total four- volume to a specific value V4 by a Gaussian distribution peaked at V4. 
In this way one avoids the presence of a zero-mode of the volume fluctuations, 
whose removal had hampered the extraction of short-distance continuum physics 
in previous work. 

After reconfirming our previous result for the de Sitter volume profile we 
measured the probability distributions V n (N) for finding a three- volume N at a 
fixed, given time step n. Well inside the universe, we found them to be Gaussian 
around the average, confirming their semiclassical nature. Moving toward the 
'stalk', where N becomes small, the probability distributions become asymmetric 
and display discretization artifacts, in the sense that the data in each V n (N) can 
be seen to separate into three distinct families. These effects are confined to very 
small length scales, but - depending on the type of observable under study - they 
can potentially interfere with the nonperturbative physical short- distance effects 
we are keen to explore, for example, along the lines suggested in [24]. 

The volume profile is an example of an observable which is not particularly 
sensitive to these discretization artifacts, even down to short scales. We found 
further evidence for this remarkable robustness when refining our spatial slicing, 
effectively doubling the number of slices in a given time interval. This was done 
by identifying the connected (on the dual lattice) layers of (3,2)-simplices in 
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between successive (4,l)-layers with slices of fixed half-integer time. Determining 
the volume profile for the (3,2)-layers alone, we found that after rescaling their 
volume by a constant p, the resulting profile could be combined with the previous 
data for the (4,l)-slices to yield a single, universal volume distribution. This 
analysis could be extended further by splitting each time interval [n, n + 1] not 
just into two but into three, yielding a similar outcome. Our results show that the 
macroscopic geometry associated with the semiclassical limit of CDT quantum 
gravity is not closely linked to the microscopic piecewise linear geometry of the 
building blocks, in line with the nonperturbative character of this limit. Our 
attempt to learn more about the nature of the A-C phase transition by studying 
the scaling behaviour of the constant p as a function of the bare parameters k 
and A yielded at this stage inconclusive results. 

Our new method of volume-fixing also enabled us to make new precision mea- 
surements of the effective action for the three-volume fluctuations 8N(n) of the 
universe around its semiclassical limit. Specifically, we investigated correction 
terms to the action - obtained by inverting the covariance matrix (8N(n)8N(n')) 
- which are of higher order in powers of A" 2 / 3 (n) and are associated with 
curvature-squared terms in the continuum. We observed nonvanishing correc- 
tion terms of this type in the discretized action, while others, expected from 
a comparison with the continuum, gave little or no contribution. This neither 
proves nor disproves the presence of higher-order curvature terms, since we found 
that their coefficients in the discrete effective action must have a specific, non- 
trivial scaling behaviour (as a function of the UV cut-off a) in order to survive 
in the limit a — > 0. We undertook a careful, quantitative analysis of one of these 
coefficients, but did not find evidence for the required scaling, at least not for 
the range of coupling constants considered. The implication is that there is cur- 
rently no evidence for the presence of higher-order curvature in CDT quantum 
gravity, although this issue needs to be re-examined once we have improved our 
algorithms to penetrate to shorter distance scales by performing the simulations 
closer to the putative UV fixed points in the phase diagram. 
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Appendix 



In this appendix, we repeat the analysis of Sec. 5 for the system where we have 
subdivided all time steps into two, as already described in Sec. 4. That is, we 
associate constant half- (odd-) integer times to the layers made up of (3, 2)- and 
(2, 3)-simplices. We want to investigate whether this finer-grained system can 
also be described by a discretized semiclassical effective action. If so, the action 
we have already derived, based on the data from integer-valued constant-time 
slices alone should be obtainable from the new effective action by integrating out 
the (3,2)- and (2, 3)- "degrees of freedom". In this way the system described by 
the integer-time slicing may be understood as arising from a "Kadanoff blocking" 
in time of the larger, finer-grained system. 

After measuring the covariance matrix of "three- volume" fluctuations 13 for all 
integer and half-integer times, we invert this matrix to obtain the effective action, 
up to quadratic terms in the fluctuations. The large inverted covariance matrix 
can be decomposed into blocks describing the (4,l)-system (at integer time), 
the (3,2) + (2,3)-system (at half-integer times) and off-diagonal blocks describing 
interactions between the two. As expected, we observe a constant shift by 2e in 
the (4,l)-block only, which we subtract from the entries of the fluctuation matrix. 

The functional form of the matrix entries resembles that found in Sec. 5, but 
with some modifications. Before proceeding further, let us rescale the number of 
(3,2) + (2,3)-vertices with the factor p determined in Sec. 4 to obtain the number 
N(n+ 1/2) which is part of the universal volume distribution N(n), valid for both 
integer and half-integer n. The measured matrix structure is consistent with an 
effective action of the form 



13 We put "three-volume" in quotation marks because - unlike the number of (4,l)-simplices 
- the number of (3,2)- and (2,3)-simplices in a slab between times n and n + 1 does not have a 
direct interpretation in terms of a three-volume. 




n 



±- (f(N(n), N(n + 1)) - ^(JV(n + 1/2)) - V 2 (N(n)j) 



n 



+ e 
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where 



and 



= i£^( 1 + Xl (^) +Xt ('-Z»)\... ] , (30) 

x + y \ \x + yj \x + y; 
x + y V \x + yj 



V 1 (N{n + 1/2)) = fi^N(n + 1/2) 1 / 3 - \^N(n + 1/2), (31) 
V 2 (JV(n)) = -^N{nf^ + A ( JiV(n). 

We observe that the "kinetic" , finite-difference terms couple not only neighbour- 
ing layers (in this case n and n + 1/2), but also directly the (4,l)-layers at n to 
those at nil, in the latter case with a negative sign. We have checked that 
this feature persists when we perform a further "time refinement" by associating 
the (3,2)- and (2,3)-simplices with distinct spatial slices, as was done in Sec. 4 
above. There are further unexpected minus-signs in the potential V 2 , and a new 
correction term linear in (x — y) in the kinetic term g(x,y). We take this as 
an indication that - unlike what happened in our earlier analysis of the volume 
distribution itself - we run into short- distance lattice artifacts for measurements 
of the volume fluctuations when considering the subdivided slicings. 

Demanding that integrating over the (3,2)-fluctuations should give back the 
(4, Ineffective action leads to explicit relations among T, I\, T 2 , as well as among 
u, and \ e g, A^, A^. A simple, but tedious computation yields 

hwrk (32) 

and 

» = (-/^ (1) + p- 1 V 2 f ^r 2 , (33) 

\ - ( \W 1 n-M 2 ^ 2r * r2 r 

as well as relations between Xi, Xi an d 

We have checked numerically that starting out with the larger covariance 
matrix of the subdivided system and integrating out the (3,2)-fluctuations re- 
produces within measuring accuracy the covariance matrix determined earlier by 
considering integer times only. The results at the point (/to, A) = (2.2, 0.6) in 
phase C are presented in the table below, which compares the directly measured 
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Figure 8: Dependence of the parameters r\ (lower curves) and r 2 (upper curves) 
on hq (left figure) and A (right figure) for selected points on the phase diagram. 



parameter 


direct 
(from (4,1)) 


integrated 
(from (4,1) and (3,2)) 


r 


23 ± 1 


26 ±2 




13.9 ±0.7 


11.3 ±0.7 




0.027 ±0.003 


0.027 ±0.003 



parameters of the (4,l)-action (15) with their counterparts obtained from mea- 
suring the fine-grained system and using formulae (32) and (33). In the fits we 
have neglect higher-order corrections, and the errors are estimated as statistical 
errors of the fits. In view of this the agreement between the two sets of values 
appears quite satisfactory. 

We have also examined the dependence of the action parameters Tj on the 
bare coupling constants, using the same values as in the analogous investigation 
of the parameter p in Sec. 4. For constant A we observe a marked dependence 
on k (see Fig. 8), in the sense that both Ti and T 2 seem to diverge when ap- 
proaching the critical line for the A-C phase transition. For fixed kq we observe a 
nontrivial behaviour of r 2 , which may be corroborating evidence that the semi- 
classical scaling function (6) ceases to give a good description of the three- volume 
distribution around A w 0.3. 
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